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thermochemical mass transfer number 
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characteristic exhaust velocity 

skin friction coefficient in presence of 
blowing 

skin friction coefficient with no blowing 
Stanton number 

specific heat of nonvaporizing components 
in solid grain 

specific heat 

port diameter 

thrust 

mass flux of gas per unit area = m /A 

g P 

mass flux per unit area at X = 0 
enthalpy 

total enthalpy (including chemical enthalpy) 

heat of gasification of solid grain 

binder enthalpy of vaporization 

thermal conductivity of gas 

mass fraction of nonvaporizing components 
in solid grain 

fuel mass fraction of surface 

oxidizer mass fraction at boundary layer 
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mass flow rate 

mass flux from solid grain 

number density of particles in gas 
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- temperature deep in grain 

- effective radiation temperature 

- surface temperature 

- X component of velocity 

- y component of velocity 

- axial distance from head-end 

- distance normal to grain 

- radiation path length 

- empirical radiation parameter 

- boundary layer thickness 

- sensible enthalpy difference between flame 

and wall 

- total enthalpy difference between flame and 

wall (including chemical enthalpy) 

-(^s - "o) 

- emissivity 

- dimensionless flame height 

- viscosity 

- oxidizer-to-fuel ratio, 0/F 

- density 

- bulk density of vaporizing component 

- variable — property correction 



0 - Stefan-Boltzmann constant 
Y - shear stress 

- velocity ratio, u/u^ 

Q - stoichiometry parameter, 

e 



Subscripts 

b - conditions at flame 
e - conditions at boundary layer edge 
f - solid grain 
g - gas 

r - in the presence of radiation 
w - wall (surface) 
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ABSTRACT 



This report presents a summary of hybrid 
rocket internal ballistics. Discussions are 
primarily concerned with the current theory 
that is based on a heat transfer limited model. 
Areas of application, major limitations, kinetic 
effects, and controversial aspects are briefly 
discussed. Areas are identified that require 
additional analytical and/or experimental in- 
vestigation . 

The review was completed in June 1971, and 
does not include bibliographic material after 
that date. 



vii 



FORWARD 



The intent of this report is : (a) to present 

a bibliography of pertinent literature related to 
hybrid rocket internal ballistics, (b) to indicate 
which references are the most significant with 
regard to the various areas of hybrid combustion 
(heat transfer limited model, kinetic effects, de- 
sign practices, etc.), (c) to present a brief out- 
line of the current internal ballistics theory, 

(d) to discuss the areas of application, the major 
limitations and controversial aspects of the cur- 
rent theory, and (e) to indicate those areas which 
require additional experimental and/or analytical 
investigation. Literature on hybrid rocket internal 
ballistics which had appeared through June 1971 has 
been included in this report. 

This work was sponsored by the Naval Air 
Systems Command, assistant commander for Research 
and Technology, Propulsion Administrator, AIR-330 
under job order number 4505-1519, 045011-01. 
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1.0 INTRODUCTION 



Although there are various types of "hybrid rockets" discussed in the 
literature, the term normally infers a system consisting of a solid fuel and 
a gaseous and/or liquid oxidizer. So-called reverse hybrids (solid oxidizer, 
liquid or gaseous fuel), tribrids , etc., have also received significant 
attention. 

Interest in hybrid rocket propulsion has fluctuated considerably since 
its inception some thirty-five years ago. However, most of the investiga- 
tions which have contributed significantly to the current understanding of 
the internal ballistics of hybrid rockets have been conducted within the 
last ten years. The majority of the contributions to this field have come 
from France (ONERA) and the United States. 

The hybrid rocket offers some significant advantages. Among these 
advantages are: deliverable specific impulse as high as that for solid or 
liquid rockets, high density impulse, space storability, ease of thrust modula- 
tion and engine restart, normally free of combustion instability, the ability 
to provide good performance with smokeless exhaust, and inherent safety. It 
is the latter characteristic, inherent safety, that has recently revived 
interest in hybrids. In addition, air-augmented hybrids show promise for 
future tactical applications. 

Many aspects of hybrid technology are similar to those of solid and 
liquid rocket engines: liquid injectors, tankage, valves, positive expul- 

sion devices, solid grain casting methods, prediction of theoretical per- 
formance, and etc. 

The one major area of difference is that of internal ballistics. 

Solid Propellant burning rate is normally characterized by an empirical 

• XI 

equation of the form r = ap and is in many cases sensitive only to pressure 
and initial grain temperature. A complete understanding of the combustion 
mechanism(s) of solid propellants is lacking and there is still considerable 
controversy in this regard. In contrast, the regression rate of the solid 
grain in a hybrid rocket may be a function of the pressure, the temperature, 
the grain configuration, and the total mass flux of gas in the port per unit 
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of cross-sectional area. Thus, predicting hybrid internal ballistics is 
dependent on an understanding of both combustion and flow phenomena. 

A bibliography of hybrid literature was published in 1967 (1)*. In 
addition, reference Cl lists much of the early experimental work (thru 1962). 
There also are several good publications which have summarized the state of 
the art in hybrid internal ballistics and in hybrid rockets in general. 
Within the last three years several additional publications have appeared 
which have a significant bearing on the current understanding of combustion 
phenomena. Although hybrid internal ballistics theory has progressed to the 
stage where it has been used for the successful design (an subsequent use) 
of several practical hybrid motors, there remain many significant problem 
areas which require further investigation. As the utilization of hybrid 
rockets increases, they will inevitably be required to perform in regimes 
of operation in which the current internal ballistics theory will be of 
limited utility. 



^Numbers in parenthesis refer to references. 
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2.0 GENERAL DISCUSSION 



A complete description of the internal ballistics of a combusting 
system would include: (a) a model for the combustion process; (b) a model 

for the internal flow (droplet history, boundary layer growth, pressure loss, 
etc.) which, when coupled with the combustion model, would yield the burning 
rate (regression rate for solid grains) ; (c) a description of the combustor 
(grain geometry, L*, contraction ratio, etc.) and how it would change with 
time; and (d) a means of predicting the ideal performance. It is the first 
two aspects that distinguish one system from another. 

The internal ballistics of hybrid rocket motors are somewhat more com- 
plex than those of solid rocket motors. This results not from a lesser 
understanding of the hybrid system but rather from the fact that the regres- 
sion rate of the hybrid grain is, in part, dependent upon the total mass 
flux per unit area passing through the port. This mass flux, in turn, varies 
with axial position due to the consumption of the solid grain. 

An understanding of hybrid internal ballistics is dependent upon an 
understanding of the interrelationships between the momentum and energy 
transport and the chemical kinetic mechanisms. A generalized sketch of the 
processes which may exist in hybrid combustion is shown in Figure 1. 

The bulk of hybrid internal ballistics modeling has dealt primarily 
with the transport of heat by convection and radiation from a diffusion 
flame zone (located in a turbulent boundary layer) to the solid surface. 

This transport is complicated by the presence of high rates of wall mass 
addition. In most of these investigations the wall gasification processes 
are considered to respond instantaneously to the heat flux and the kinetic 
processes are neglected. Recent investigations have modified this viewpoint. 
Other investigators have studied the solid grain regression rate from the 
alternate approach of the thermal degradation processes of the solid. 

Neither approach is complete in itself. The former has been more popular 
because it relates directly to the physical design variables and resulted 
in a successful engineering model which is applicable to many current hybrid 
propellant systems. However, future propellants and applications may render 
the current heat transfer limited model inadequate. 
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It will be very difficult to handle simultaneous heat and mass transfer 
together with kinetic mechanisms in a working model for hybrid internal bal- 
listics design. Fortunately, in many current and proposed hybrid systems it 
appears that the "limiting case" models will be applicable. 

Hybrid internal ballistics modeling has yielded good qualitative re- 
sults and in limiting cases has yielded quite accurate quantitative results. 
Currently, any quantitative inaccuracies or regimes of operation outside the 
model assumptions are handled through empirical correction parameters. 

Hybrid internal ballistics literature is voluminous and of varied con- 
tent. Some of it emphasizes the development of analytical models whereas 
others emphasize the engineering applications. Most of the early (pre 1965) 
modeling of hybrid combustion has been reviewed by Green (54) . Comprehensive 
development of the current internal ballistics model resulted primarily from 
investigations at the United Technology Center (UTC) . However, other inves- 
tigators at the Lockheed Propulsion Company (LPC) , Rocketdyne, Stanford 
Research Institute (SRI), and etc. have made pertinent contributions. Some 
of the more important references are listed here by corporate authors. 

LPC: C13 (2 volumes), C14, C15 

MIT: 38, 39 

ONERA (France): 20, 22, 25, 26 
Rocketdyne: 30, 31, 32 

SRI: 17 

UTC: 3, 4, 5, 7, 8, 10, 12, C4, C5, C6, CIO 
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3.0 BRIEF REVIEW OF HEAT TRANSFER LIMITED MODELS 



3 . 1 Introduction 

The material presented in this section is primarily that developed by 
UTC although reference is also made to other material when it is appreciably 
different. No attempt has been made to present all of the details and justi- 
fications for the model development; these can be found in the cited refer- 
ences. However, some detail is presented to indicate major aspects so that 
subsequent discussion can be made with regard to the applicability and limi- 
tations of the model. 

In brief, the model developed in this section has proven to be appli- 
cable (although with some reservation) to many of the current hybrid rocket 
propellant systems and regimes of operation. Thus, the hybrids that can be 
described reasonably well by the model are those which have fairly simple 
grain configurations (no sharp corners, no flow mixers within the grain, 
etc.), which are pure hybrids (little if any solid oxidizer in solid grain), 
and which operate at relatively high pressures. 

The basic model considers the solid phase regression rate to be con- 
trolled and limited by convective and radiative heat transfer to the solid 
grain from a relatively thin diffusion flame in a turbulent boundary layer. 
The existence of this physical description has been verified by schlieren 
observations (13) and studies with boundary layer probes (8, 9). All chemi- 
cal reactions are considered to occur very rapidly in an infinitely thin 
flame sheet and therefore, all kinetic effects (thermal degradation of 
binder, heterogenous attack of solid by active oxidizer at the wall, gas 
phase kinetics, et.) are neglected. The possible existence of a liquid 
metallic melt on the surface is also neglected in the following development 
but will be discussed subsequently. The basic model is applicable only to 
the steady-state and considers the surface to respond to the heat flux 
instantaneously by decomposing to yield gases (and possibly solids). Thus, 
in the steady-state 

6 = m h = p^rh = (pv) h [1] 

w f V f V wall V 
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where 



h = 

V 

r = 

i = 

Pr = 



heat of gasification of the solid grain 
regression rate of solid grain 
total heat flux to the surface 
density of the solid grain 



3.2 Convective Heat Transfer (2, 3, 5, 7, 8, and 12) 

The convective heat transfer to the wall can be written 



where 



% = t2] 

^ w 



Q = convective heat flux to surface 
c 

= average specific heat of gas 
k = thermal conductivity of gas 
h = enthalpy 

y = distance normal to grain surface 

A Stanton number is defined for heat transfer from the flame to the wall in 
the presence of mass addition such that 

Q. 



where 



Lh = fc - fc T ) 

V P /b V P /w 

Ah = sensible enthalpy difference between the flame and the wall 
u^ = velocity at the flame 
= gas density at the flame 
[3] into [1] with Q,, ~ yields 



[3] 



[4] 



= h 



Q 0, 

— = — = C p u, Ah/h 
h H b b V 

V V 



[5] 
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An expression is now needed for Cj^ in the presence of blowing. Empirical 
expressions for the skin friction coefficient, C^, in the presence of blowing 
appear in the literature and it is therefore convenient to determine 

Cjj = f(cp [6] 



Assuming that the Reynolds analogy is applicable between the flame and the 
wall in the presence of blowing (or Le = 1) yields for constant but nonunity 
Prandtl number (8, 53) 



where 




9u 
P 9y 



[7] 



h^ = total enthalpy (including chemical enthalpy) 



T = shear stress 



Pr = Prandtl number 



Le = Lewis number 
Near the wall (2 , 8) 




[ 8 ] 



[7] into [8] and integrating between the wall and the flame yields 



Qw 



= Jlpr--67 



For no combustion between the flame and the wall Ah^ = Ah. 
[9] into [5] with Ah^ = Ah yields 



[9] 



T Pr-« I 

Q = I flame to wall for heat transfer [lOJ 

H 2 I 

^b^ ) 



The skin friction coefficient is defined in the conventional manner by 



T 

W 




C.P 

fee 



Boundary layer edge to wall for 
momentum transfer 



[11] 
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[11] into [10] yields 






2 

p 

e e 
pb^b 



Pr 



-.67 



[ 12 ] 



[12] into [5] yields 



Pfr 






p u B Pr 
e e 



where 



B 






a thermochemical parameter for a 
particular propellant combination 



[14] 



It remains to find an expression for in a turbulent boundary layer with 
mass addition at the wall. 



For no blowing (66, p. 537) 




.03 Re 



- 0.2 



[15] 



Marxman (5) has developed an expression for which can be approximated 

by the following expression over the indicated range of the blowing 
parameter B ' . 

— « 1.2(B' (5 < B' < 100) [16] 

So ~ 



where 



B’ 






(17J 



Smoot, et al (C14) used Spalding's thin film model for the effect of blowing 
on wall friction. In that case 



C 



f 



C 



fo 




[16a] 



A comparison between Marxman 's exact solution (5) and the approximate 
expressions [16] and [16a] is presented in reference C14 and indicates 
that [16a] is better than [16] for B’ <5 and that for B* from 10 to 25 
they are equivalent. Most hybrids have been found to have a value of B’ 
between 7 and 15. 
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[15] and [16] into [13], 



p^r = .036 



G Re Pr 

X 



[ 18 ] 



where 



G = ^0^0 ” total mass flux of gas per unit area 



[19] 



and 



Re = 



GX 



e e 



[ 20 ] 



However, from [13] and [17] with p .r = (pv) 

f w 



B' = B Pr 

Thus with Pr = 1 [18] can be written 



-0.67 



[ 21 ] 



• - 0.2 0 23 

p.r = .036 G Re B 

f X 



.0.36G»-8(^r‘’'\»-23 



[ 22 ] 

[23] 



Equations [22] and [23] apply to incompressible flow and neglect variable 
fluid properties in the boundary layer. Marxman (5) has accounted for 
variable fluid properties in a semi-empirical manner and obtains 



P 



f 



r = 




[24] 



where 

p is a reference density and (p/p^)^*^ generally has a value between 
1.0 and 1.1 but is difficult to determine accurately. 



To this point the following major assumptions have been made: 

(a) flat plate turbulent boundary layer analysis is applicable. 

(b) regression rate of grain is controlled by heat transfer from a 

diffusion flame. 

(c) flame zone is infinitely thin. 

(d) oxidizer enters port as a uniform gaseous stream. 

(e) no heat transfer into subsurface region of solid grain. 
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(f) Reynolds analogy is applicable between the flame and the wall in 

the presence of blowing. 

(g) Le = Pr = 1 

(h) Pressure is constant. 

(i) no radiation heat transfer to the wall. 

It has been analytically determined (8) that the axial pressure gradient, 
which results from the mass addition (with no particle lag) induced flow 
acceleration, increases by approximately 10%. It was also shown (8) 
that particle lag losses (metalized grains, etc.) and the pressure gradient 
have only a negligible effect on the shear stress distribution and the 
velocity profile. Thus, the flat plate analysis is applicable with little 
error . 

The blocking effect which reduces normally depends only upon the 
gas mass addition at the wall since solids which leave the surface usually 
occupy negligible volume compared to the gas. If the hybrid grain has solids 
which do not vaporize, equation [1] should be written as follows (8) 



Qw 



— = p 



r 



f V 



P h 

V V 



eff 




+ C AT 
m 




[25] 



where 

T^ - T^ = surface temperature - temperature deep in the grain 
specific heat of nonvaporizing component 

^v _ mass of vaporizing component _ bulk density of 
V volume of grain vaporizing component 

_ mass of solid additives 
m^ total mass of grain 




Thus , the binder is heated from temperature T^ and vaporized and the solid 

additives are heated from T to T . Equation [24] should now be written as 

o s 

follows 



Q 

P^r = p^(l - K)r = — 



eff 




[26] 
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where 



B 



u 

e 



Ah 






h 



eff 



[27] 



For hybrid systems which satisfy the above assumptions, equation [26] 
provides a working equation for internal ballistics calculations. However, 
it remains to be shown that B, given by equation [27], can be evaluated. 
Various combustion models have been employed in order to evaluate B. The 
various combustion models (5, 43, 47, 54, etc.) have been discussed by 
Marksman and Wooldrige (4). Marksman (5) assumed Le = Pr = 1 (similar 
profiles of enthalpy, concentration, and velocity) and used a Prandtl 
mixing length approach to solve for B, Q , (}) , and 



where 



n = 



Kf V 

s 

K 

ox 



V = 






K 



ox 

e 

n, = 



ox-fuel mass ratio 

fuel mass fraction at surface 

oxidizer mass fraction at boundary layer edge 

velocity ratio = — 
u 

e 



= dimensionless flame height 



The complete solution (5) reveals that for all practical purposes B does not 
vary with axial distance, X (assuming that L/D < 25 so that severe oxidizer 
diluation does not occur). The mixture ration, v, and the oxidizer mass 
fraction at the boundary layer edge, K , were shown to decrease with in- 

OXg 

creasing X. The concept of a nons toichiometric diffusion flame has caused 
considerable debate (see discussion by Green in reference 4) and points 
out the present lack of understanding of transport processes in gaseous 
combustion. 
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3.3 Radiative Heat Transfer (3, 5, 7, and 8) 



Radiation to the surface of the hybrid grain is treated as a grey-body 
equilibrium problem in which the flame zone (i.e. , the gas or dense particle 
cloud) is considered as a radiating continuum. It has been pointed out (7) 
that this simplified approximation, although presently found to be adequate, 
may be inadequate in future systems which use more advanced propellant sys- 
tems with high radiation heat transfer. However, improvements in this area 
will require a better understanding of nonequilibrium radiation than is 
presently available. The radiant heat flux is written 



Q = oe 
r 



wO 



r- 

r 




Now let 



where 



^ oe e 

w g r 




= effective radiation temperature 
a = empirical constant 



[28] 

[29] 

[30] 



z = radiation path length 



N = p for gas systems, or 

= n, the number density of particles, if there are radiating 
solids in the gas 



Then , 

Reference 17 presents a derivation of an expression for e^ wiiich is con- 
siderably more complex than equation [30]. Although the derivation yields 

a better understanding of the parameters affecting e , the resulting expres- 

§ 

sion does not appear to offer any distinct advantage over equation [30]. 



A -aNz\ 

(1 - e ; 



(31) 



3.4 Combined Convective and Radiative Heat Transfer 



The regression rates of most metalized hybrid grains will result from 

• • 

both convective and radiative heat transfer. The effects of and can- 
not be summed directly since is a function of the surface mass addition. 
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Increased heat flux to the wall increases r which increases the surface 



blowing. This, in turn, decreases Q but has little effect upon Q . 

/ \q 5 ^ ^ 

Neglecting fp/p 1 * and letting Pr = 1, equations [18] and [31] yield. 



= h 






= .036 G Re 



n o / \ n 77 (-1 -aNz'\ 

“0-2 /„,\-0.77 „ w r U - e ) 



X 



eff 



(b;) 



B + 



[ 32 ] 



eff 



The thermochemical parameter B is not significantly effected by Q but the 

* 

blowing parameter B’ will increase with increasing Q^. 



(pv) 



B’ = B’ with radiation = 



w 



P ^ 

V 



p u C /2 p u C /2 
e e fr e e fr 



[33] 



[33] into [32] yields 



h B’p u /2 = h .036 G Re 
V r e e fr v __ x 

eff eff 



-0.2 (b;^ B + 



[34] 



From [13] 



[without radiation = ^C^/2^ ^e^e^ ^ ^ 

eff 



-0.2 ^0.23 ^ 

D n 

X V 



[35] 



Dividing [34] by [35] and noting that 



fr 






eff 



[36] 



yields 



B’ 

C ^ 






[37] 



[37] may be approximated by (3) 



B’ 

r 



1.3 *^r 



[38] 



From [32] 






, \ - 0*77 



eff 



eff 



[39] 



14 



[38] into [39] yields 



V =r 



-Q /Q 

r c , 'r 
e + 



eff 




[40] 



where 6 is the convective heat transfer without radiation, 
c 

Equation [40] indicates the coupling that occurs between the convective 
and radiative heat transfer. Two systems may have the same regression rate 
but if is significantly present in one, then it will respond differently 
to changes in pressure, oxidizer mass flux, and etc. From equation [31] it 
is observed that the radiation heat transfer is not directly dependent upon 

G as is the convective heat transfer. However, increased G (which increases 

• • 

r) does affect by increasing the pressure and/or the number density of 
radiating particles (C6) . Reference C14 (which used the Spalding thin film 

model) arrives at a somewhat different result than equation [37]. The 

• • 

coupling is handled by a graphic technique (without making an approxi- 
mation corresponding to equation [38]) and the results indicated that equa- 

• • • 

tion [40] underestimates r for coupled systems with B < 100. However, 

equation [40] presently appears to be adequate and more readily lends itself 
to an internal ballistics model which can be programmed for engineering design. 



In summary 



V - Pf ■ h 



V 



-Q /Q 

r c 



eff 



[41] 



Q = .036 h 

c V 



eff 




0.6 




- 0.2 



g0.8b0.23 



[42] 



G = 




[43] 

[44] 



These equations, when coupled with grain geometry and theoretical performance 
data provide the necessary material for an internal ballistics computer pro- 
gram which can be used for engineering design and analysis. 
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3.5 Parameter Evaluation 



In order to use equations [41] through [44] for internal ballistics 
calculations several parameters must first be evaluated. Inspection of the 
equations indicates that the following parameters must be evaluated: h. 



^Vb’ ^s’ ^P/Pe) » n, T^, and z. 



eff ’ 



Details of the approximations and/or calculations used in order to ob- 
tain these parameters can be found in the literature (8, 12, C4 , etc.) and 
will only be briefly considered. 



h - the effective heat of gasification of the solid phase. 

^ef f 

This parameter must be calculated or determined experi- 
mentally. In order to determine h analytically, 

veff 

equation [25] indicates that the surface temperature, , 
the effective heat of vaporization of the binder, h , and 

vb 

the specific heat of the nonvaporizing solids, C^, must be 
evaluated. The latter parameter is normally known or can 
be estimated from the literature. 



h^ - the effective heat of vaporization of the binder. The 

^ thermal decomposition processes of the binder must be 

understood if h is to be analytically determined. 

^b 

Calculation of h^ for simple polymers such as plexiglass 

have been made (12, etc.). For more complex binder systems 

or solid grains which incorporate solid oxidizers, h^ will 

normally have to be evaluated experimentally, h^^ may be 

significantly reduced due to surface reactions caused by 

the oxidizer in the grain. In addition, for those systems 

in which oxidative attack of the binder is important 

(discussed below) , one method of accounting for this 

pressure sensitive effect would be through a pressure 

sensitive h (or h ) . 

^b ^eff 



T - the surface temperature. T is normally estimated to be 
s s 

between 600 and 800°K (39). 
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- the effective radiation temperature. This parameter is 
estimated with any error absorbed in the empirical radia- 
tion parameter a. UTC (8, 12, etc.) has used = 2/3 

(T ^ , X. « ) • radiating systems, the effect of 

stoichiometric ® 

the local mixture ratio, 0/F, can be significant. As noted 

above, O/F is found to decrease with increasing distance 

from the head-end. Any initial dilution of the oxidizer 

(N 2 in air, etc.) will also significantly reduce T^. A 

reduction of 0/F also occurs when a hybrid motor is 

throttled and this will affect T . These effects can be 

r 

readly included in the current model by treating T^ as 
T^ (x,t) or T^ (0/F). 



B - the mass transfer number. Equation [27] indicates that 

evaluation of B requires evaluation of u /u, , Ah, h 

e ’ veff 

As pointed out in references 8, 12, etc., even significant 

errors in the determination of Ah and h do not sig- 

^eff 

nificantly effect the calculated r in convective heat 

• 0.23 

transfer limited systems since in those systems r B 
In order to find (5, 17) it is necessary to evaluate 

(p/p^) , which may be determined analytically (5). In UTC’s 
current internal ballistics program (12), .036 (p/p)^*^ in 
equation [42] is treated as an empirical constant which is 
determined experimentally . 



0/F - the local mixture ratio at the flame. The value of 0/F 

at the flame is known to be fuel rich (3, 32). UTC (8, 12, 

etc.) estimates 0/F to be (3/4) (0/F) ... . and 

stoichiometric 

considers it to be independent of X. However, as noted 
above, 0/F may vary with axial position, oxidizer, dilution, 
and throttling. Although this does not normally affect 
significantly it may affect (i.e., T^) as indicated 
above. It appears that a variable 0/F could be incorporated 
into the current model without too much difficulty. 
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Ah - the sensible enthalpy difference between the flame and 

wall. Once the local mixture ratio and the surface tem- 
perature are known or estimated, Ah can be readily 
calculated (10), 

- the effective wall emissivity, is estimated. Muzzy 

(12) has used e = 0.90. 

w 

a - the empirical radiation parameter, a is determined from 
experimental data. 

- the effective port area (that yields an average u^) . 
is the actual port area less the area occupied by the 
boundary layer displacement thickness (12), 

m - the total mass flux of gas in the port, m varies with 
g g 

the axial position, X, It must account for the mass flux 

of gas entering the port , the mass flux from the decompos- 
ing solid grain, and any gas consumed in the production or 
comsumption of solid particles (10), In addition, it should 
be emphasized that m^ includes only the vaporized oxidizer 
(CIO), The model assumes uniformly injected vapor. How- 
ever, in real systems the time required for complete vapor- 
ization of the oxidizer may be significant. This correc- 
tion could be included in the current model but depends 
upon knowledge of the initial drop size distribution which 
is not readily determined. 

n - the number density of particles in the flame zone. The 
contribution of the radiation from the gas is neglected 
when radiating particles exist in the flame zone. Expres- 
sions have been developed for calculating n (8, etc.) 
depending upon the source of the particles. The resulting 

expressions show that n varies as p/m and is a function 

§ 

of the axial position, X. In the evaluation of n, the 
gas temperature (T^) at the boundary layer edge, must be 
determined. An approximate method for evaluating T^ has 
been presented (8, 10, 12, etc,). 
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z - the optical path length. This parameter is not well under- 
stood since actual systems have nonequilibrium radiation. 
Reference 8 has considered z proportional to the hydraulic 
diameter of the port for radiating gaseous systems and as 
a constant (proportional to flame height) for radiating 
particle systems. Any inaccuracies in n or z are accounted 
for in the empirical radiation parameter a. 

Pr, 

, etc. - these parameters are calculated by normal techniques, or 
estimated . 

As indicated above, the regression rate model must be coupled with 

grain geometry equations and performance parameters such as the ideal thrust 

coefficient, C , the ideal characteristic exhaust velocity, C* , and the com- 
r 

bustion and nozzle efficiencies. These parameters are evaluated in the same 
manner as for liquid or solid propellant systems. Several grain configura- 
tions and propellant combinations have been evaluated. References C4 and 
C5 present calculated and/or measured values of B, Ah, etc., for various 
propellants. In addition, the literature contains many references which 
present theoretical performance calculations for various hybrid propellant 
combinations (22, etc.). 

3 . 6 Typical Theoretical Results of the Coupled Hybrid Model 

Figure 2, from reference 12, indicates the type of results which are 

• • 

obtained from the coupled model for a typical aluminized grain. At 

low values of G , radiation heat transfer dominates the total heat flux to 
o 

the surface and therefore, the regression rate becomes practically inde- 
pendent of mass flux through the port but very sensitive to changes in 
operating pressure and mixture ratio. At high values of , convective heat 

transfer dominates and the regression rate becomes insensitive to pressure 

0 8 

changes but varies as G ’ . When the mass flux is at some intermediate 

• • 

value, both and are significant and the regression rate depends upon 
pressure, mass flux, port diameter, and grain length. The slight increase 
in r at low predicted for radiating systems as shown in Figure 2 results 
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from the increase in n (n « l/m ), and therefore Q , at low G (10). However, 

g r o 

the model assumes is independent of and, if is allowed to vary as 

discussed above, the form of the curve at low G may change. Although some 

experimental data exist (40, 41) for Q at low G , much additional experi- 

r o 

mental and analytical work is required. 

3. 7 Other Heat Transfer Limited Models 

As mentioned above, the heat transfer limited hybrid internal ballistics 
model has been developed into a practical design tool by both the United 
Technology Center (UTC) and the Lockheed Propulsion Company (LPC) . The dis- 
cussion above dealt primarily with the model as used by UTC. The LPC work 
(C13, C14 , C15) is very similar. Although LPC handled coupling, 

C^ = C^ (B*), and condensed surface species in a somewhat different manner, 
practical application of either treatment remains essentially the same. 

Other investigators have used somewhat more simplified analyses for 
heat transfer limited hybrid operation (50, 51, etc.). Recently Drew (52) 
has used the basic UTC model for systems in which a metallic melt layer is 
formed on the surface of the grain. In the analysis it is assumed that the 
melt layer does not affect the gas-phase heat transfer. Lieberherr (26) has 
also considered the effects of a metallic melt layer on the surface. 

Various investigators from Japan have also studied hybrid combustion 
(59 , 74) . Reference 59 used a Rayleigh flow model in place of the boundary 
layer treatments discussed above. These authors treat the flow as laminar 
which is normally not applicable to hybrid environments. 
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Fig. 2 REGRESSION RATE DEPENDENCE ON OXIDIZER MASS 
FLUX AND PRESSURE FOR A TYPICAL ALUMINIZED 
HYBRID SYSTEM (from p.51 Ref. 12) 
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4.0 KINETIC EFFECTS IN HYBRID COMBUSTION 



4. 1 Introduction 

As mention above, the heat transfer limited model has been found to be 
adequate for many propellant systems over restricted mass flux and pressure 
regimes. However, experimental data at low pressure and high mass flux (8, 
17, 27, 28, 29, etc) indicate that the heat transfer limited model is not 
applicable in these regimes for many propellant systems. Smoot and Price 
(29) have found that at high G, low P, the regression rate becomes independent 
of G but is directly related to the partial pressure of the oxidizer in the 
gas phase. In addition they have found experimentally (27, C14) that the 
addition of metal to the solid grain reduced the pressure dependence of r. 

There has been considerable debate concerning which kinetic mechanism(s) 
in the overall combustion process become rate controlling in these limiting 
cases. The various mechanisms will be discussed below but first it is im- 
portant to briefly consider whether or not kinetic effects need to be 
studied from the practical standpoint of expected operational environments. 

It is true that for a wide range of propellant systems that the regres- 
sion rate is dominated by heat transfer over much of the expected range of 
mass flux. In these regimes of operation the heat transfer limited model 
can be expected to be adequate for internal rates and minimum grain surface 
area the trend in some future hybrid systems will undoubtably be towards 
the highest mass flux possible. The practical upper limit on G^ or G de- 
pends on ignition and flooding effects which will be discussed below. 

However, at these high mass fluxes, kinetic effects may become significant. 
The propellant system could still be characterized by sub-scale testing 
and/or the heat transfer model could be modified to include an empirical 
correction term which is pressure dependent (10) , but neither of these solu- 
tions provide any understanding of the basic phenomenon. This understanding 
will be required if hybrid systems of the near future are to be adequately 
modeled and scaled. 

There are several other developments which require an understanding of 
the kinetic effects in hybrid combustion. Solid oxidizers are being in- 
cluded in the fuel grain in order to obtain higher regression rates and 
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smaller variation in I with 0/F. As the percent of oxidizer is increased, 

sp ^ ’ 

the hybrid system will become more like a conventional composite solid pro- 
pellant in which kinetic mechanisms dominate, A somewhat common limiting 
problem exists for hybrid systems with significant oxidizer in the solid 
grain, as well as for solid propellants subjected to high gas velocities 
parallel to the surface - i.e., in erosive burning environments. An under- 
standing of the interrelationships between the kinetic mechanisms and the 
turbulent transport mechanisms would be of great value in both of these 
closely related areas. 



There is currently a considerable interest in air-breathing hybrids 
and solid-fuel ramjets. The heat transfer model discussed above can be 
expected to provide the basic model for these systems since oxidizer dilu- 
tion has been considered in some detail for hybrids (assuming the grains 
are not kinetically dominated and if the boundary layer flow is not signifi- 
cantly altered) . When the oxidizer is diluted the flame moves further away 
from the surface, and increases in thickness (3), while the flame temperature 
drops (C4). An analysis by UTC (C4) has shown that the primary effect of 
dilution in Q dominated systems is on the parameter AH (and therefore B) . 
Dilution may also significantly effect r in a dominated system through 
the parameter T^. 

Significant kinetic and radiation effects can be expected in many air- 
breathing applications and these aspects have received only limited attention. 

A recent analysis (discussed below) has indicated that surface catalytic 
attack by active oxidizing species at the grain surface may be the control- 
ling mechanism in the high G, pressure dependent regime of hybrid operation. 

If this proves to be true, then it may be found that as the oxidizer is sig- 
nificantly diluted, pressure sensitive combustion may occur at significantly 
lower mass fluxes than found in conventional hybrids. Variable thrust hy- 
brids, in which the oxidizer is throttled, may also present an operational 
environment in which the combustion mechanism is controlled by heat transfer 
at high thrust (high pressure) and controlled by kinetics at low thrust 
(low pressure) . 

High frequency combustion instability has been observed only rarely in 
conventional hybrids. Combustion instability will be discussed below but is 
mentioned here because of its obvious link to kinetically controlled combustion. 
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When kinetic mechanisms contribute significantly to the hybrid combustion 
process (air-augmented systems, high oxidizer content in solid grain, low 
pressure from throttling, etc.) combustion instability may become a sig- 
nificant problem. 

It is apparent that the kinetic aspects of hybrid combustion may be of 
considerable importance in the near future and an understanding of these 
phenomena is required. As mentioned above," there has been considerable 
debate in the literature concerning whether the kinetically influenced 
regime of combustion occurs in the gas phase at the flame or heterogeneously 
at the wall. A brief review of the literature related to the kinetic aspects 
of hybrid combustion is presented below. 

4 . 2 Kinetic Modeling 

Until recently most investigators have considered that the pressure 
sensitive regime of hybrid combustion was due to either gas phase kinetics 
or exothermic heterogeneous reactions at the surface. In addition, several 
investigators have studied hybrid regression rates by considering solid 
phase kinetics. 

Smoot and Price (27, 28, 29, C13) found experimentally that the regres- 
sion rate of LiH/rubber fuels was a function of and the partial pressure 
of active oxidizer but practically independent of total pressure. Their 
experimental data also indicated that burning rate catalysts had no effect 
on pressure sensitivity and that the addition of metal reduced the pressure 
sensitivity. They reasoned from this data that the pressure sensitivity of 
r was controlled by the turbulent transport of active oxidizer to the wall 
and exothermic first order heterogeneous reaction of this active oxidizer 
with the solid fuel. Solid fuel decomposition, radiative heat transfer, 
and binary gas phase reactions were not considered to play a dominant role. 

An analytical development is presented in reference C13. Houser and Peck 
(31) have also indicated that oxidative degradation is required to explain 
the high regression rates of hybrid fuels. In contrast Miller (42) con- 
sidered the pressure sensitivity to result from reduced reaction rates at 
low pressures which broaden the gas phase flame zone. 
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Kosdon and Williams (33) have attempted to explain r = r(p) by consider- 
ing small perturbations to the heat-transfer limited regime. Marxman ( 8 , 

17) has discussed the limitations of this model. An interesting approach 
has been presented by Tsuge and Fujiwara (59) who have argued that a Rayleigh 
flow model better fits the hybrid situation of high rates of wall mass addi- 
tion than the boundary layer model. Their treatment was restricted to laminar 
flow and kinetic behavior was introduced through an Arrhenius expression for 
surface pyrolysis. The model correctly indicated the independence of r on 
G at high values of G. 

Rabinovitch (48) studied the thermal degradation of polymers. His study 
indicated that the wall temperatures found in hybrid rockets are adequate to 
produce the observed hybrid regression rates through thermal degradation 
alone. Recently Lengelle (40) has also studied thermal degradation kinetics 
but the results were inconclusive. Marksman ( 8 , 8 , 17) has also analyzed 
pressure sensitive hybrid combustion and, supported by Rab inovitch ’ s work, 
has argued that gas phase kinetics in the flame zone are the principal con- 
trolling mechanism during low pressure, high G^ operation. His development 
(7) indicated that the pressure sensitive regression rate cannot be explained 
solely by heterogeneous surface kinetics. Subsequently, Marxman, et al (17) 
have conducted experiments in the pressure sensitive regime. Two binders 
(PU, PBAN) , three oxidizers (O 2 , F 2 > N 2 O 4 ) , and various aluminum loadings 
(0 , 20, 40%) were employed. Tests were conducted with varying G^ and varying 
throat area. Schlieren studies also showed that the flame zone broadens at 
low pressures. The experimental results indicated that the percent aluminum 
and the binder decomposition did not significantly affect the pressure depen- 
dence of r (although considerable data scatter existed at the lowest pres- 
sures) . The former result for aluminum effects is in contrast to the results 
of Smoot and Price for LiH (28). Lieberherr’s (26) work suggests that the 
pressure sensitivity of hybrid grains containing metallic hydrides may be 
due in part to surface melting and liquid flow. 

Marksman, et al (17) reasoned that "at high pressures the gas-phase 
reaction rates are fast enough to consume all of the available material that 
is being vaporized at the surface." It is further reasoned that gas-phase 
reaction rates decrease with decreasing pressure until the flame cannot 
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consume all vaporizing material. This results in a shift to a mixture ratio 
which is even more fuel rich, and thus is reduced. The reduction in 
decreases r. An analytical expression was derived for the pressure sensi- 
tivity of the steady-state regression rate. The turbulent diffusion flame 
which exists at high pressures was considered to broaden and transition to 
a turbulent pre-mixed flame at low pressures. By using wrinkled flame 
theory (64) for the premixed turbulent flame- and Denison and Baum's (63) 
result for the kinetic flame speed, an expression is developed for r in the 
pressure sensitive regime. The resulting expression adequately fit the ex- 

perimental data for the pressure sensitivity of r at a fixed G but did not 

• ^ 

properly predict the independence of r on G at high G. Additional experi- 
mental and analytical work is needed to determine the validity of this 
model . 

Recently Stickler and Kumar (38, 39) have studied the thermal degrada- 
tion of polymers. Primary consideration was for polymethyl methacrylate 
(PMM) . They have shown that thermal degradation cannot by itself account 
for the high regression rates of hybrids with wall temperatures within the 
range reported in the literature (600-800°K). It is argued that gas phase 
kinetics are not slow enough and that exothermic heterogeneous reactions 
require more oxidizer at the wall than reported experimentally. An analysis 
of turbulent transport through the flame zone is presented and the conclusion 
is reached that the small percent of active oxidizer that diffuses through 
the flame to the wall acts as a catalyst for depolymerization. Thermal 
degradation reduces the polymer to large fragments and the catalytic action 
is considered to break down the fragments to a size which has a vapor pres- 
sure at the surface equal to the chamber pressure. The rate of catalytic 
depolymerization is considered proportional to the oxidizer concentration 
at the wall. It is argued that at high pressures the catalytic rate is 
sufficiently high that the burning process becomes heat-transfer (thermal 
degradation) limited. At lower pressures and high G, the regression rate 
is limited by oxidizer transport to the wall. In this limit there is not 
enough oxidizer at the wall to adequately degrade the polymer for vaporiza- 
tion and r becomes practically independent of G. The authors have pointed 
out that the kinetics of oxidative degradation are not well understood and 
therefore, it will be difficult to actually predict regression rates in the 
pressure sensitive regime. 
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The work of Stickler and Kumar appears to have merit although the 
observed broadening of the flame zone at low pressures (17) , the lack of 
effect of various binders on the pressure sensitivity (17, 27), and the 
possible pressure sensitivity of the formation and consumption of metallic 
melt layers (26) indicate that pressure sensitive hybrid combustion is a 
complex process which probably involves more than one pressure sensitive 
mechanism. 

The above discussion indicates that a considerable lack of understanding 
still exists in the area of pressure sensitive hybrid combustion. Although 
several new analytical models have appeared in the past few years, very few 
basic experimental studies have been conducted. Additional experimental 
investigations are required. The inherent data scatter obtained at low 
pressures needs to be considered further. A real need exists for an experi- 
mental method that allows accurate determination of the regression rate 
during combustion. 

In addition, the similarity between the treatments for the rates of 
thermal degradation of polymers (38, 30, 48, 40, etc.) and the rates of 
fracture of polymers subjected to compressive or tensile loads (67, 68, 69, 
70, etc.) requires investigation. 
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5.0 APPLICATION, IMPLICATIONS, AND LIMITATIONS OF 
THE HYBRID INTERNAL BALLISTICS MODEL 



5 . 1 General Discussion 

Although the current internal ballistics model does not include kinetic 
effects, it can be used for many propellant combinations over a fairly wide 
range of mass flux. There are many important design questions which can be 
answered by utilizing the heat transfer limited model. In addition, there 
are areas of hybrid operation which are outside the basic model assumptions. 

It is not the purpose of this section to discuss in detail all of the 
applications and limitations of the current model which have been considered, 
but rather to indicate some of the major considerations and to reference the 
literature where pertinent discussions have been presented. The discussion 
below is arbitrarily divided into three sections: Ignition and Mass Flux 

Considerations, Regression Rate Tailoring and Injector Effects, and Per- 
formance. These topics are obviously interrelated. 

5 . 2 Ignition and Mass Flux Considerations 

With the natural trend toward the use of high oxidizer mass flux, it 
becomes necessary to better understand the internal ballistics in this 
’’limiting" regime of operation. Each propellant combination has an upper 
and lower mass flux beyond which operation becomes irregular and the bal- 
listics model is not applicable. 

As discussed above, for high mass fluxes many propellant systems become 
dominated by kinetic effects. However, even those systems that remain in- 
sensitive to pressure into the high mass flux range can reach flow rates in 
which ’’Flooding’’ occurs. Ignition becomes a problem at high and the 
burning process may exhibit periodic pressure spikes or extinguish completely. 
Wooldridge and Muzzy (10) have discussed come of the problems related to 
ignition and operation at high mass fluxes. In the case of hypergolic pro- 
pellant combinations, high G^ operation can cause excess oxidizer to accumu- 
late over a portion of the grain which may subsequently ignite and cause 
severe pressure spikes. In non-hypergolic propellant applications the high 
G may prevent adequate oxidizer vaporization after ignition and combustion 
may be extinguished. 
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Various methods of hybrid ignition have been considered by Muzzy, et al 
(CIO), Hypergolic propellant combinations, solid propellant pyrotechnics, 
solid propellant ignition boosters to enhance ignition (C9) (which degrade 
the safety advantage of hybrids) , and hypergolic liquid bipropellants have 
been used. 

Ignition lags from several hundred milliseconds down to tens of micro- 
seconds have been observed (C7, 22). Various techniques have been employed 
to decrease the ignition lag (C8, etc.). In general, it is advantageous to 
minimize ignition lags in order to obtain higher reliability and repeati- 
bility and to prevent ignition pressure spikes (22) . Ignition overpressure 
is highly undersirable. In addition to the degradation or destruction of 
structural integrity, these overpressures can cause an initially high mixture 
ratio (0/F) which may result in significant throat erosion in the first few 
seconds of operation (C9). This throat erosion in turn can significantly 
degrade the overall performance and decreases the utility of the internal 
ballistics computer program. 

Very few studies of ignition phenomena in hybrids have appeared in the 
literature. Barrere (22, 25) has initiated work in this area and other in- 
vestigations are also pertinent (34, 75, C12 , etc.). The ignition process, 
both initially and in restart operations, is probably one of the least under- 
stood phenomena in hybrid combustion. Many design variables appear to effect 
the ignition process. Currently many questions remain partially or fully 
unanswered. For example, what is the best way to distribute the ignition 
energy and/or the oxidizer flow over the hybrid grain (uniformly, head end, 
in valley, on plateaus, etc.)? The ignition energy can also be distributed 
in various forms — gases, solid-liquid interaction, etc. What is the optimum 

amount of ignition energy consistent with the propellant formulation, G , 

o X 

etc? Can the hybrid tolerate overignition? What actually causes the 
ignition? What effect does initial oxidizer temperature have on ignition 
(C4) ? What is the optimum time between ignition and oxidizer flow initia- 
tion? Will addition of solid oxidizer to the solid grain insure greater 
ignition reliability and less ignition lag? These are some of the problems 
which must be considered. There is an obvious need for basic studies in 
ignition which would hopefully lead to at least a qualitative model for the 
hybrid ignition process. 
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At the other extreme of mass flux-low G -a '’limiting" G also exists, 

o ° o 

but for considerably different reasons than at high G^. Problems associated 
with low mass flux and the resulting low regression rates of the solid grain 
appear when variable thrust operation is considered. Boost-sustain duty 
cycles or continuous throttle capability often necessitate very low G , low 
r operation. In addition, when restart capability is required, heat soak 
into the solid grain between firings must be considered. Ignition at low 
G^ may also be expected to produce undersirably long ignition delays. 

At very low G^ several model assumptions may be violated. The turbulent 
boundary layer may transition to laminar and this may explain the significant 
data scatter typically found at low G^ (C14) . In addition, heat may penetrate 
to significant depths below the surface resulting in a nonsteady regression 
rate. 



Significant pressure effects may also be present at low G in metalized 

o 

systems due to the pressure sensitivity of the radiative heat transfer. 

The problems associated with transient operation (ignition, restart, 
throttling) and operation at low G^ , low r have received some analytical 
attention by several investigators (3, 7, 8, 10, 50, etc.). After ignition 
or throttling, a period of time is required to establish the steady-state 
burning rate (temperature profile in the solid grain). This time may vary 
in restart systems depending on the soak period between firings (50). It 
has been shown (8, 10) that this time is critical for solid grains which 
have a high vaporization temperature or a low enthalpy of vaporization. 
Operation at low r and/or in the transient regime may cause subsurface 
material to be chemically or physically altered (charred, melted, decomposed) 
which in turn can effect both steady-state combustion and restart capability 
(C16) . Operation in this regime may also result in low frequency pressure 
oscillations associated with alternate accumulation and removal of the melted 
or charred layer on the surface (17) . 

Current practice is to run subscale tests to determine the upper and 
lower mass flux limits for a particular propellant combination. Incorporat- 
ing transient phenomena directly into the current internal ballistics pro- 
gram would be difficult. Subscale data and results from simplified transient 
analyses can be used as "corrections" to the ballistics model during the 
transient portion of the duty cycle. 
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The internal ballistics of throttled hybrids within the heat transfer 
limited operating regime requires additional considerations beyond those for 
throttled liquid systems. This occurs because of the coupling between r and 
G or G^. A brief consideration of the implications of this coupling is pre- 
sented below. Details can be found in references 10, C4 , C6 , CIO, and etc. 

Various methods of obtaining thrust-time (F - t) variations in hybrids 
have been discussed in the literature (10, C4 , CIO, etc.). Programmed F - t 
missions can normally be handled through appropriate grain design; geometry, 
multiple castings, etc. When demand F - t is required, some form of oxidizer 
throttling is normally required. Throttling of head-end oxidizer is the 
simplest method (provided adequate performance and stability can be obtained 
without aft injection) but has a distinct disadvantage. When G^ is reduced 
G^ decreases, but only to a much lesser extent since r varies approximately 
as Thus, as the system is throttled the pressure drops, the 0/F de- 

creases, and in many cases this would result in decreased performance. 

For demand thrust, reference C4 has discussed several alternatives: 

(1) provide excess propellant so that the motor can operate off optimum O/F, 
and sacrifice the performance; (2) use fore and aft-end oxidizer injection in 
which the combustion products would be fuel rich for all head-end G^. Aft- 

end G would then be varied to obtain the desired thrust at the optimum 

o 

mixture ratio (10) ; (3) use a propellant system with a pressure sensitive 
regression rate so that the fuel flow rate varies more directly with the 
oxidizer flow during throttling. However, operation in this pressure sensi- 
tive regime may introduce other difficulties such as combustion instability, 
etc.; (4) develop a propellant with a large value of n where n is the expo- 
nent in the regression rate expression r = High values of n would 

allow throttling at practically constant 0/F but present problems with re- 
gard to matching grain progressivity with the resulting rapid decrease in 
• 

r that would occur as the grain was consumed. Multiple port designs with 
separate oxidizer injectors have also been used (CIO) . Throttling is accom- 
plished by shutting off some of the injectors. 

The throttling capability of hybrids is one of its inherent advantages 
over conventional solids. It is difficult to estimate what throttling ratio 
can be achieved without undue loss of performance or complexity of design. 
This will depend on propellant choice and future techniques for throttling. 
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5 . 3 Regression Rate Tailoring and Injector Effects 



Because of the interdependence of most of the operating parameters in 
hybrid internal ballistics it is difficult to separate discussion of the 
various aspects. Two important areas that are isolated in this section are 
regression rate tailoring and the problems from injector effects, but the 
discussions are keyed to interrelationships with maxx flux and ignition. 



5.3.1 Regression Rate Tailoring . One of the problems with hybrid propul- 
sion systems is the very low regression rate which characterizes most pro- 
pellant combinations. Various techniques are needed for increasing the rate 
(r) and for obtaining time variations. Good detailed discussions can be 
found in references 8, 12, C4 , C5 , and C15. 



From equations [27] and [42] it is observed that for convective heat 
transfer limited systems 




[45] 



For heat transfer limited systems in which both convective and radiative 
effects are present the situation is more complex and equations [41] through 
[44] must be considered simultaneously. In their applicable ranges these 
equations indicate what can and cannot be done to effect a significant in- 
crease in f. In addition, they point to other techniques (with reservation) 

which may be employed that are beyond the limitations of the model. 

• 

In Q - limited systems, nominal changes in Ah and h^ rr have little 
c ^ / \23 eft 

effect upon regression rate since p^r « f^^/^^veff )* ’ Therefore, addition 

of exothermic material to the solid grain to reduce h^ cannot be expected 

, eft 

to significantly increase r (C4) . However, by adding solid materials that 
do not vaporize at the surface, but rather leave the surface as solids, 
is decreased and r increased correspondingly. There is an obvious practical 
limit to the amount of solids which can be added. If adequate overall per- 
formance is to be obtained solids must be efficiently burned in the gas 
phase (8). It should also be noted that for the upper limit where all solids 
are consumed in the gas phase, a significant amount of oxidizer may be con- 
sumed above the flame, thereby reducing the partial pressure of the oxidizer. 
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In systems which are dominated by radiative heat transfer, r may be 

significantly effected by changes in Ah or h^ since surface blowing does 

err ^ 

not counter the heat transfer as readily as it does in systems. 

The above discussion has indicated several means by which the regres- 
sion rate can be increased; addition of highly exothermic material to the 
solid grain, addition of nonvaporizing solids to the solid grain, and in- 
creasing the radiative contribution to the overall heat trnasfer to the 
grain. The latter is normally accomplished by increasing the metal content 
of the grain which increases the effective radiation temperature, T^, The 
maximum metal content is limited by the amount that can be efficiently con- 
sumed in the overall combustion process. In addition, high metal loading 
may cause a melt layer to form on the surface. As mentioned above, the 
metallic melt layer may be periodically blown off, giving rise to low fre- 
quency pressure oscillations or undesirable combustion characteristics (C16) , 
The formation of the metal melt also requires some modification to the bal- 
listics model as discussed by Drew (52) and Lieberherr (26) , 

Another method of increasing r is by the addition of solid oxidizer 
to the grain. The quantity must be limited to prevent self-sustaining com- 
bustion if hybrid operation is to be maintained. Reference C4 has indicated 
that this limiting quantity depends both on the type and size of the solid 
oxidizer crystals. Addition of solid oxidizer will also introduce kinetic 
effects and the regression rate will become pressure dependent with increased 
possibility of combustion instability. The gases leaving the surface and 
entering the diffusion flame probably consist of binder and oxidizer decom- 
position products and the products from fuel rich combustion near the surface. 
Consequently, the quantity of solid oxidizer that renders the current heat 
transfer ballistics model inadequate is likely to be considerably less than 
the maximum allowed for non-self sustaining operation. 

Reference C4 has also discussed the possibility of adding materials 
to the grain which yield gas-producing reactions below the surface. This 
would cause mechanical ejection of solid materials into the gas stream where 
they must then be consumed. 

It has been found that port configurations that induce mixing produce 
higher average regression rates. Reference C4 reports that a triangular 
port produced a significant increase in r over the conventional circular port. 
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The addition of an orifice (diaphragm, restrictor) to the grain to induce 
large-scale turbulent mixing has been found to increase the regression rate 
downstream of the orifice by as much as 50% (17, 22, 24). Various additives 
have been found which significantly increase r at high pressures and produce 
pressure sensitive combustion (15) . 

It is interesting that the heat transfer limited model can point to 
methods by which r can be increased but at the same time it is not applicable 
in many of the resulting combustion situations. 

5.3.2 Injector Effects . Injector effects are one of the more important 
aspects of hybrid combustor design. Injector and head-end design can sig- 
nificantly affect the oxidizer vaporization rate and oxidizer droplet pene- 

• 0 .8 

tration of the boundary layer. Since r « G * , injectors that gouge out the 
solid grain at the head end (due to wall impingement of liquid or vaporized 
oxidizer or to severe recirculation zones) may significantly alter the re- 
gression rate throughout the entire grain and alter the expected performance. 
These effects can easily overshadow all ballistic model calculations. Large 
injector effects have been observed for lengths up to five port diameters 
from the injector. In these regions the boundary layer is disturbed and the 
ballistic model is not applicable. 

There is a real need to study injector effects and to model this 
process. This will be difficult and perhaps empirical correlations are the 
most that can be expected. Ideally, it would be desirable to be able to 
account for oxidizer break-up, oxidizer droplet vaporization and/or com- 
bustion rate, droplet trajectories, and interaction of the oxidizer drops 
with the surface of the grain. Except for the latter, these have been 
attempted in liquid rocket studies with only limited success. In fact, the 
inability to accurately predict droplet size distributions produced by in- 
jectors remains to be one of the major problems which prevent meaningful 
quantitative predictions with vaporization-limited combustion instability 
models. In boost phase operation at high pressures, some liquid oxidizers 
will be in environments exceeding critical conditions near the head-end of 
the grain and this can significantly change the required analysis for the 
consumption rate of the oxidizer drops. The interaction of liquid oxidizers 
with combusting or decomposing solid surfaces is also not well understood. 
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The injector effects are directly related to the limiting opera- 
tions discussed in Section 5.2. Injectors which do not permit adequate dis- 
tribution and vaporization of the oxidizer may cause local surface flooding 
with the resulting phenomena discussed above, and in addition, irregular 
consumption of the solid grain and low performance. 

Limited attempts have been made at analytically predicting oxidizer 
droplet histories in hybrids. LPC has briefly considered the effects of 
drop size and velocity, and motor geometry on oxidizer vaporization (C13). 

Because of the current lack of adequate modeling of injector effects 
in hybrids, most solutions to injector related problems have been solved 
experimentally by trial and error. UTC (C4) has used splashb locks of a low 
regression rate fuel at the head end to prevent severe head-end erosion, to 
enhance vaporization, and to allow the vaporized oxidizer to enter parallel 
to the walls of the grain. Uniform regression rates have also been obtained 
by inserting a diaphragm (orifice, restrictor) near the head end of the 
grain (17, 22,24). 

5.4 Performance 

All of the above discussions are of course directly related to obtain- 
ing performance from the hybrid rocket. The limits of upper and lower oxi- 
dizer mass fluxes and the 0/F shift that occurs as a result of throttling 
were discussed above and must be considered when attempting to obtain maximum 
performance. In this section several additional considerations are discussed. 

Because r a X design of high performing hybrid grain con- 

figurations is more complex than for solids. As the grain is consumed G 
decreases (for a mixed injected mass flow rate) and therefore r decreases. 

If constant thrust and high performance are to be maintained, a progressive 
burning area will be required that matches the reduction in r. There is 

competition between two desirable hybrid features (10) . In order to minimize 

• 

r changes with time, changes in G need to be minimized. This necessitates 
a large initial port size. However, it is also usually desirable to have 
the maximum volumetric loading possible which in turn means that the port 
size should be small. A compromise must be reached which depends upon the 
specific application (10) . 
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To minimize performance losses which result from 0/F shifts (due to 

throttling or grains which are not progressive enough) , propellants with 

flat I vs 0/F curves should be employed. Reference C6 indicates that in 
sp 

general, low cost propellant systems characteristically have a sharp peak on 

the I vs 0/F curve at (I ) , and high density propellants character- 

sp ma^ 

istically exhibit the desired flat vs 0/F. It has also been reported 

that the addition of small amounts of AP to the solid grain tends to flatten 

the I vs 0/F curve (C6) . 
sp 

Because some of the vaporized fuel is convected downstream below the 
flame, it has been found necessary to provide mixing of the gaseous products 
in order to obtain high performance (C4, C5 , C13, 24, etc.). Various tech- 
niques have been considered and employed. UTC (C4) (and subsequently other 
investigators) added gas stream turbulators or mixers aft of the grain in 
an aft plenum region and obtained significant increases in delivered I 
ONERA (24, etc.) (and subsequently other investigators; (17, etc.)) used a 
diaphragm (orfice, restrictor) within the grain port and various boundary 
layer tripers have also been used. The diaphragms and boundary layer trips 
evidently destroy the turbulent diffusion flame downstream of the device by 
creating high turbulence which in turn has yielded high regression rates. 

The heat transfer limited boundary layer model is of little use for internal 
ballistics calculations in these applications. The well-stirred reactor 
concept (44, 45, 46, etc.) is in this case probably a better model. 

Aft-end oxidizer injection provides mixing and high performance. 
Multiport grains provide high fuel mass fluxes and provide good mixing at 
the aft end of the grain. Submerged nozzles have also been used to aid 
mixing (8) . 

Reference C13 has shown that although low performance results pri- 
marily from poor mixing of the gases, liquid metal on the surface and poor 
evaporation of the liquid oxidizer also contributed to performance losses. 
Lieberherr (26) has also discussed the effects of a metallic surface melt 
on performance. 

Aside from the important basic considerations of grain progressivity . 
mixing of the gases, and maintainence of a constant mixture ratio are several 
design aspects which may overshadow the utility of the ballistic model (s). 
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Injector effects, which were discussed above, are one such design aspect and 
can result in poor fuel utilization, low performance, and unpredictable 
thrust-time behavior. Throat erosion can also be significant (especially 
during high pressure, high operation) and can result in significant de- 
gradation in performance. The rate of throat erosion is difficult to pre- 
dict but in many cases it must be establish for incorporation into the bal- 
listics program if meaningful design calculations are to be obtained. 

It has been shown (3, 8) that in many cases r does not vary appreciably 

with distance from the head-end. For L/D > 5 the boundary layer edge becomes 

the port centerline. Oxidizer dilution then occurs but as long as it is not 

severe (L/D < 25), B has been shown to remain practically constant (5, 8, 

etc.). The increase in G with X tends to offset the decrease in r with X 
f * 08 23 —0 2^ 

(r a G * B* X * Thus, for the convective heat transfer limited 

hybrids, r does not vary appreciably with X. If the model is to be used 
for L/D > 25 then B must be considered a variable. In highly radiative 
systems the shift in 0/F due to oxidizer dilution may significantly affect r. 



Because of the numerous variables available to the hybrid designer, 

past experience will often be required in order to make final decisions on 

alternate designs. Excellent performance (% theoretical C*, % theoretical 

I , fuel utilization) from hybrid propulsion systems has been reported by 
sp 

numerous investigators (22 , C5 , C7 , C8, C13, etc.). 
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6.0 COMBUSTION INSTABILITY 



Studies of combustion instability in hybrid rockets have been pri- 
marily experimental (4, 16, 17, 22, C4 , etc.). Occurrences of combustion 
instability in hybrids appears to be rather infrequent to date and, when it 
has occurred, was normally related to injector failures, inadequate vapor- 
ization, etc. Occurrences of acoustic type coupling have been observed only 
very rarely. However, with solid oxidizers being added to the grain, opera- 
tion in low pressure regimes or with large amounts of oxidizer dilution being 
contemplated (ducted rockets, etc.), and with turbulators being used to pro- 
mote a well-stirred reactor within the port, it should be expected that the 
occurrence of combustion instability will increase in the future. Combustion 
instabilities can result in increased regression rates, system overpressuri- 
zation, and depending on the type of instability, may either increase or 
decrease the combustion efficiency. 

Reference 17 has discussed three general types of combustion in~ 
stability that may occur in hybrid rockets. The first type of instability 
is related to the vaporization and combustion of the oxidizer drops and is 
normally of low frequency. This type of instability is somewhat related to 
that found in vaporization-limited liquid propellant rockets but is probably 
more complex. Liquid rocket combustion instability is still not well under- 
stood and much debate continues with regard to the mechanism(s) by which the 
instability is driven (71, 72). Some models for liquid rocket combustion 
instability are concerned with the coupling between the vaporization process, 
a pressure disturbance, and the fixed combustor geometry. Instabilities of 
this type are normally high frequency acoustic modes. In hybrid rockets, 
motor geometry changes with time. The liquid droplets vaporize similar to 
that in liquid rockets (although probably much slower) but in most cases the 
vapor is not as rapidly consumed. The vapors must be transported by turbu- 
lent diffusion to the fuel vapor which is evolving from the grain surface. 
Liquid oxidizer may also penetrate the boundary layer and interact with the 
grain surface. This may cause periodic localized flooding of the grain with 
the resulting pressure spikes discussed above. Elimination of this type of 
instability has been accomplished through injector modifications which change 
the spray and energy distribution in the port. Reference C4 reports that 
elimination of the fuel rich regions in the aft plenum eliminates pressure 
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fluctuations. Multiport designs were successful in this respect. Moutet 
and Barrere (22, 23, discussion in 4) have eliminated this type of insta- 
bility by using a diaphragm in the port as discussed above. This device 
increases the residence time of the droplets and restricts them to the head- 
end. Although this device appears to suppress the low frequency instabilities 
which result from oxidizer droplet vaporization lags, reference 17 reports 
that it had no suppressing effect upon acoustic modes of instability. 

The second type of instability is also of low frequency and many 
times non-periodic. It is "associated with periodic accumulation and break 
off of char layers or melted layers (Possibly metallic) at the surface" (17) . 
This phenomenon has been discussed above with reference to low - low r 
operation. This type of instability has been reported by various investi- 
gators (C4, C5 , C13, etc.). Elimination of this type of instability may 
require one or more of the following; reduction in metal content, change 
of binder type, use of aft-mixers, use of nonsymmetic ports, and etc. (4). 

The third type of instability is related to solid propellant combus- 
tion instability phenomena and could be expected to occur most frequently 
in hybrids which have kinetically generated pressure sensitive regression 
rates. Kinetic aspects of hybrid combustion (Section 4.0) are therefore 
of major importance. There is considerable debate with regards to combus- 
tion instability mechanisms in solid propellants (73, etc.) and acoustic 
mode instability in hybrids will probably be even more difficult to charac- 
terize and model. Reference 17 has reported the occurrence of spontaneous 
instabilities in the longituTiinal mode and also possibly in a transverse 
mode. From the discussion in Section 4.0 it is obvious that many factors 
may influence acoustic mode instability: oxidizer dilution, solid oxidizer 

in grain, pressure level, port configuration, mass flux through the port, 
and etc. 

The model developed by Marksman, et al (17) for the steady-state 
regression rate pressure sensitivity was discussed above (Section 4.0). 

That investigation was primarily concerned with modeling acoustic mode 
combustion instability in hybrid rockets. Their experiments indicated 
among other things that (a) the metal content of the solid grain did not 
affect the pressure sensitivity of the regression rate or the frequency of 
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the instability, (b) binder pyrolysis characteristics did not affect the 
regression rate pressure sensitivity, and (c) use of more reactive oxidizer 
resulted in less pressure sentitivity. 

Only limited experimental and analytical (7, 17, 18, 37) investiga- 
tions have been conducted in this area and much additional study is required. 

The use of hybrid rockets for tactical weapon systems may necessitate 
studying acceleration effects for both the liquid oxidizer droplet behavior 
and the regression rate of the solid grain. Both of these areas have re- 
ceived only brief consideration (Cll, C13) . 
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7.0 SUMMARY OF AREAS OF APPLICATION OF 
‘ HEAT TRANSFER LIMITED INTERNAL BALLISTICS MODEL 



The internal ballistic model developed by United Technology Center 
is found to be applicable — 

• With very good accuracy to propellant systems which are 
convective heat transfer limited and which have relatively 
simple grain geometries. Minimum subscale testing is 
required . 

• # 

• With good accuracy to propellant systems which are heat 
transfer limited but with both convective and radiative 
contributions to the total heat flux. Subscale testing 
is required to evaluate empirical constants. 

In both cases the models are limited to hybrid operation between the 
upper and lower critical mass fluxes. Scaling of hybrid combustors has been 
successfully demonstrated (12, C5 , C6 , CIO). Throttling applications are 
readily handled unless transient phenomena are of long duration. The major 
items which overshadow the calculations made with the ballistics model are 
injector effects and throat erosion. 

In addition, the model can be used — 

• With lesser accuracy and scaling reliability for systems 
which have slightly pressure sensitive regression rates 
as a result of kinetic effects. Additional empirical 
parameter required in model. Additional subscale testing 
is required. 

• Probably with decreased accuracy for some air-augmented 
hybrids and solid fuel ramjets if care is used in de- 
fining the mass-flux employed in the model and if regres- 
sion rates are not kinetically controlled. 
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8.0 SUMMARY OF REQUIRED AREAS OF INVESTIGATION 



1. Basic investigations are needed in the area of turbulent trans- 
port raechanisms in boundary layers with combustion and high rates of wall 
mass addition. 

2. Experimental and analytical investigations of combustion insta- 
bility in hybrids are required. These are required in order to determine 
whether or not combustion instability will actually be a future problem in 
hybrid combustion. Studies of the regression rate and oxidizer vaporiza- 
tion rate responses to pressure and velocity disturbances are required. 

3. Directly related to combustion instability studies are required 
studies of the kinetic aspects of hybrid combustion. This remains a 
controversial subject and will have increased importance as more advanced 
hybrid systems (AP added to grain, etc.) are considered. 

4. Hybrid ignition is a current problem that requires immediate 
attention . 

5. Related to ignition phenomena are studies required on injector 
ballistics effects. This is also a current problem. These studies will 
probably have to be empirical in nature. 

6. Radiation heat transfer in hybrids (and in air-augmented systems) 
requires further study. The current simplified approach used in the bal- 
listics model may not suffice for higher energy propellant systems which 
may have significantly increased radiation transport. 

7. Investigations are required to more adequately determine the 
applicability of the current model to air-augmented hybrids and solid fuel 
ramjets . 

8. Improvements can be made in the current ballistics model 
(variable mixture ratio and radiation temperature, unvaporized oxidizer 
effects , etc . ) . 
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